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Abstract 

We study the masses of the Higgs bosons in the most general two-Higgs-doublet 
model in a basis-independent approach. We adapt the recently developed Minkowski- 
space formalism to this problem and calculate traces of any power of the mass-matrix in 
a compact and reparametrization-invariant form. Our results can be used to gain insight 
into the dynamics of the scalar sector of the general 2HDM. 

1 Introduction 
1.1 General 2HDM 

The two-Higgs-doublet model (2HDM) is one of the simplest extensions of the Higgs mech- 
anism of the electroweak symmetry breaking beyond the Standard-Model, [HOE]. In this 
model one introduces two doublets of Higgs fields, <pi and <p2, which interact with the matter 
fields and also self-interact via an appropriate Higgs potential. 

Higgs potential of 2HDM contains many free parameters, which are not constrained by 
experiment. Extensive studies conducted over past decades have shown that playing even 
with a small subset of these free parameters one can get a rich spectrum of models with 
different phenomenologies (see examples in [31 H] ) . 

Recently, it has become clear that not all of these free-parameters are equally important. 
One has a certain freedom in choosing the basis in the Higgs field space when writing the 
lagrangian. This basis change shifts the values of the parameters of the lagrangian, but by 
construction it has no effect on the physical observables. Thus, it is only the basis-invariant 
features of the theory, and not the entire set of free parameters, that really shapes the phe- 
nomenology of the model. 

When using reparametrization transformations, one is immediately led to the most general 
2HDM, whose Higgs potential contains all possible electroweak-invariant quadratic and quartic 
combinations of the two doublets. Several groups have recently focused on the properties of 
the general 2HDM and developed a set of basis-invariant tools adequate for this task. The 
motivation behind this interest is not to provide the most accurate description of the real world, 
but rather to understand the whole spectrum of possibilities offered by the second doublet. 



In this way, the most general 2HDM with no a priori restriction on its free parameters should 
be viewed as a useful tool for building specific models with predefined properties, and it is 
definitely worth studying in as much detail as possible. 

1.2 Approaches to the general 2HDM 

The main problem with the general 2HDM is that it cannot be worked out with straightforward 
algebra. The obstacle arises at the very first step: when minimizing the Higgs potential, one 
arrives at algebraic equations of high order, which cannot be solved in the general case. In 
this situation, any method that would give any non-trivial insight into the model is welcome. 

Following the early suggestion of [5] , a very elaborate basis- independent treatment of gen- 
eral 2HDM was presented in [6] and further developed in [TJ [8]. In this approach one writes 
the Higgs potential as 

V = + ^Z abcd (<Pl<Pb)(<Pl<P d ) , 

and manipulates with Y a b and Z a b c d as tensors rather than just a collection of parameters. 
Instead of finding explicitly the vector of vacuum expectation values, v a , one adds it to the set 
of objects to manipulate with, keeping in mind, however, that it satisfies the extremum condi- 
tion. Along these lines, one can find several algebraically independent invariants constructed 
as full contractions of the available tensors, and some of the properties of the model could be 
seen through the prism of these invariants. 

Unfortunately, this powerful technique lacks intuition, as the results arise from lengthy 
(and often computer-assisted) algebra of invariant polynomials. A more appealing approach 
to the general case was suggested and developed in (9j El [EE1 121 [13j [H] . In this approach 
one works not in the space of Higgs fields a , but in the real four-dimensional space of gauge- 
invariant bilinears (<fi' a <fib) (the orbit space), which has the Minkowski-space signature. Many 
of properties of the Higgs potential can be derived in a very intuitive way based on simple 
geometric considerations. 

This approach was developed further in [T5], [16], where the reparametrization group was 
extended also to non-unitary transformations of the fields. In the 1 + 3-dimensional orbit 
space it leads to the full Lorentz group of transformations, and this freedom provides even 
more insight into the properties of the 2HDM potential. In particular, many of the statements 
about the general 2HDM are much more naturally formulated in terms of the four eigenvalues 
of the Minkowski tensor A M1/ (see the next Section) rather than its space-like part. Another 
key point of [T5l [16] was that all essential results were formulated exclusively in terms of 
the parameters of the potential, without using the yet-unknown vacuum expectation values. 
For example, within this approach one could formulate conditions for the spontaneous CP- 
violation and draw the full phase diagram of the model exclusively in term of the parameters 
of the potential, without using the unknown vacuum expectation values. 

Thanks to all these approaches, we have now a fairly detailed understanding of the prop- 
erties of the Higgs potential and of the vacuum in the general 2HDM. 

1.3 Towards dynamics of the general 2HDM 

The next step in the study of the general 2HDM is to understand its dynamics. This includes 
the mass spectrum of the physical Higgs bosons, the pattern of their interactions, as well as 
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their couplings to the fermions. All this must be done within a basis-independent approach. 

Let us stress once again that if one chooses a restricted Higgs potential, for example, an 
explicitly CP-symmetric one, the entire calculation is drastically simplified. One can explicitly 
find the minimum of the potential and calculate the masses and the interaction of the Higgs 
bosons. This straightforward approach fails for the most general 2HDM, which calls upon 
more involved techniques for the analysis of its properties. 

The mass spectrum of the general 2HDM was studied in a number of recent papers. For 
example, in [Til HH EE] the mass matrix was explicitly calculated in a specific basis and not 
in a reparametrization-invariant form. An interesting study was also presented in [17], where 
certain bounds and relations between the masses and the parameters of the potential were 
observed, however that work relied only on numerical analysis. Finally, very recently a very 
detailed account of the dynamics of the general 2HDM was presented in [8]. Among other 
results, explicit expressions of the mass matrix were derived in U (2)-invariant way in terms of 
various full contractions of tensors Y a t> and Z a b c d as well as vacuum expectation values. 

In the present paper we show how to analyze the masses of the physical Higgs bosons in 
the Minkowski-space formalism. We obtain compact 5*0(1, 3) -invariant expressions for the 
traces of any power of the mass matrix. Thus, we can now use the full power of the extended 
reparametrization symmetry of the problem to gain further insight into dynamical properties 
of the general 2HDM. 

The structure of the paper is the following. In Section 2 we briefly review the Minkowski- 
space approach to the general 2HDM. The main results of the paper are then derived in 
Section 3. For each possible type of the 2HDM vacuum, we find the mass-matrix and calculate 
the traces of its powers in reparametrization-covariant way. A discussion of the results and 
conclusions are presented in Section 4. 



2 Overview of the formalism 



In this work we focus on the scalar sector of 2HDM. The Higgs potential of the most general 
renormalizable 2HDM, Vh — Vi + Va, is conventionally parametrized as 



V 2 

v 4 



"111 (01 </>l ) + "4(0202) + mi 2 (0l0 2 ) + m? 2 *(0201^ 



y (0 f i0i) 2 + y(0 2 2 ) 2 + A 3 (0 t 1 0i)(40 2 ) + A 4 (0l0 2 )(0|0O 
1 
2 



A5(0l02) 2 + A*(40l) 2 



A 6 (010l) + A 7 (02>2)1 (0 f 10 2 ) + h.C 



It contains 14 free parameters, four in the quadratic and 10 in the quartic terms, which 
makes the phenomenology of 2HDM very rich even at tree level. However, not all points 
in this 14- dimensional space of parameters lead to distinct physics: if two sets of parame- 
ters can be mapped into each other by a certain linear transformation between the doublets 
(reparametrization transformation, or basis change), they will lead to the same physics, [S1H]. 
Usually one insists that the kinetic term be invariant, so one considers only global unitary 
transformations between the two doublets, U (2). However, as shown in [UJUn], one can extend 
this reparametrization group to the general linear group GL(2, C). The Higgs kinetic terms is 
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not invariant under non-unitary transformations, but it can be treated in a reparametrization- 
covariant way, so that all the physical observables still remain invariant under this extended 
reparametrization group. This approach has provided several new insights, which would be 
very difficult to see using the more traditional unitary reparametrization group. 

Technically, the extended reparametrization group can be implemented as follows. We 
switch from the fields to bilinears and introduce the four- vector r M = (tq, r») = 
where $ = (fa, fa) T is a 2-dimensional vector of Higgs doublets and a 1 are the Pauli matrices. 
This four- vector is gauge invariant and parametrizes the gauge orbits in the space of the Higgs 
fields. The general reparametrization group GL(2, C) can be written as C* ® SL(2, C), where 
C* is the group of simultaneous multiplication of both fa with the same complex number, 
while 5*1/(2, C) is the special linear transformation group. It is the latter group that leads to 
non-trivial transformations of the Higgs potential, which we now focus on. 

Transformations of $ under SL(2,C) correspond to the SO(l,3) transformations of r M , 
equipping the gauge orbit space with the Minkowski-space structure. It follows from the 
definition of that 

r = (4<!>i) + (402) > , r% = 4 [(4fa)(4<l> 2 ) - {4cj> 2 ){4fa)} > , (2) 

so that the physically realizable vectors populate not the entire 1 + 3- dimensional Minkowski 
space, but the future lightcone (LC + ). The Higgs potential in the r^-space can be written in 
a very compact form: 

V = -M^ + ^A^rV . (3) 

Here the four- vector M M is built from parameters m?- in (TjQ), while the symmetric four-tensor 
A^;, is constructed from the quartic coefficients \. Their explicit expressions as well as some 
properties can be found in [131 EE EBL Here we just note the most important property of 
for potentials stable in a strong sensqj: A M „ can always be diagonalized by a certain SO(l,3) 
transformation of the r^-space, and after diagonalization it takes form 

\ lu = diag(A , -Ai, -A 2 , -A3) with A > 0, A > A; , i = 1, 2, 3 , (4) 

where the inequalities among the eigenvalues result from the positivity constraint on the 
potential. The minus signs in front of the "space-like" eigenvalues arise from the pseudo- 
euclidean metric in the orbits space. 

It is known that the potential ([3]) can have three types of minima: (i) the electroweak 
(EW) conserving, (ii) the EW-breaking but charge conserving (i.e. neutral), and (iii) the EW- 
and charge-breaking ones. One can use the v.e.v.s of the two doublets (fa) to construct (r M ). 
Then, the three type of minima correspond to: (i) (r^) = (the apex of the forward lightcone 
LC + ), (ii) (r>*) ^ but (r^)(rA = (the surface of LC+), (iii) (r>*) ^ and (^(r^) > 
(the interior of LC + ). The position of the charge-breaking extremum (r u ) c h is given by the 
following equations: 

A^(r u ) ch = M^ (5) 



1 We use here the terminology of [TT]: the potential is stable in a strong sense, if its quartic part increases 
along all rays starting from the origin in the Higgs field space. The potential is called stable in a weeak sense, 
if the quartic part has flat directions, but the quadratic potential increases along them. For the Minkowski- 
space analysis of potentials stable in a weak sense, see |18j . where a similar condensed-matter problem was 
considered. 
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If is not singular, a solution of this system always exists and is unique: (r^)^ = = 
(A -1 ) ^ U M U . However, the requirement that (r u ) c h lies inside the forward lightcone places 
bounds on that could yield physically realizable solutions. In addition, the charge-breaking 
extremum is minimum only if all Aj < 0, i — 1, 2, 3, i.e. if the tensor A^ is positive-definite in 
the entire space of non-zero vectors r M . When searching for the neutral extrema, we use the 
Lagrange multiplier technique. The positions of all neutral extrema (r M ) are the solutions of 
the following simultaneous equations: 

A^<7v)-C = M^, C M = C<^>, (6) 

where C is a Lagrange multiplier. This system can have up to six solutions, [TOj [11], [15] , among 
which there are at most two minima, while the other are saddle points, [16J. 



Finally, following [15], we write the Higgs kinetic term covariantly as 

K = , (f = (d a $) , (7) 

where a denotes the usual space-time coordinates, while fi, as before, refers to the orbit space. 
The reparametrization transformation properties of (f 1 are the same as for r M . In the "default" 
frame, = (1, 0, 0, 0). Upon an 5*0(1,3) transformation, acquires non-zero "space-like" 
coordinates, however the condition K^K^ = 1 is always satisfied. The four-vector is not 
involved in the search for the minimum of the potential, however it affects the mass matrix at 
this minimum. This generalized kinetic term effectively incorporates the non-diagonal kinetic 
term, which, as was argued in [T9] . must be introduced in the initial lagrangian to restore 
renormalizability of the model. 



3 Mass matrix of the most general 2HDM 

In the previous studies, [T5l[T6] the Minkowski-space formalism was used to understand various 
properties of the 2HDM lagrangian and of the vacuum state. The next logical step is to study 
the dynamics of the model in an reparametrization-covariant way. In this paper we make 
a step towards fulfilling this program. We obtain expressions for the mass matrices of the 
physical Higgs bosons in the most general 2HDM and study some of their properties. 

When doing so, we stick to the Minkowski space formalism, but we adapt it to our prob- 
lem. Although the masses are physical observables and are reparametrization-invariant, the 
mass-matrix is, obviously, basis-dependent. So, for intermediate calculations we switch back 
from the bilinears to the Higgs fields themselves, derive the mass-matrix in a specific basis, 
then calculate the traces of the powers of this matrix, and return to the Minkowski-space 
formalism. Although the resulting equations do not yield the masses in a simple closed form, 
they nevertheless can be useful for the analysis of the general 2HDM. 

In the subsection devoted to the neutral vacuum below, we also comment on relation of 
our results with some of the previous studies of the mass spectrum. 

3.1 Switching to the real fields 

Let us denote the complex fields as <p ii(X , where i = 1,2 indicates the doublet, while a =f, j 
indicates the upper and lower components in a given doublet. Let us then introduce the 



5 



8-component real vector of scalar fields tp a , a = 1, 8, with the following components: 

(p a = (Re0i iT , Im0i jT , Re0 2 ,T, W>2,t, ^ e( t>i,i^ lm 0i,h R e 02,|, W> 2 ,|) • ( 8 ) 

The four- vector r M can be rewritten in terms of <p a as 

= ipj?^ . (9) 

Here, S M are four real symmetric 8-by-8 matrices; S° is just the unit matrix, while explicit 
form of S* can be immediately reconstructed from the definitions (see Appendix). Since the 
upper and lower components of the doublets are not mixed by the Higgs potential, matrices 
S M have a block-diagonal form, composed of identical 4-by-4 matrices. Below, we will often 
deal with these 4-by-4 matrices, denoting them by the same letter S M . Which set of matrices 
is being used, 4-by-4 or 8-by-8, should be clear from the context. 

In contrast to a^, the matrices S M do not form a closed algebra, but they belong to a larger 
algebra (S^,rP), described in the Appendix. They also share with a^ 1 an important property: 

{£*,!?'} = 28 ij -Is, (10) 

where brackets denote the anticommutator. It follows then that if a regular real symmetric 
8-by-8 matrix A is written as a^'E fl , then its inverse is 

A -i = a^_ S^ = (S°, -S l ). (11) 

Below we will encounter products of matrices E's and E's. When simplifying these products, 
the following results prove useful: 

^(WZ U +E"E' i ) = cT -Is, (12) 

With this notation, we can give a compact expression for the mass matrix in a specific basis. 
Let us write the expansion of the scalar lagrangian near an extremum as 

1 d 2 V 

L ~ (K p J: p ab )(d a ip a )(d a if b ) - H ab (ip a - {<pa))(ipb - (<Pb)) , H ab = - , 

where the hessian H ab is calculated at the extremum. The 8-by-8 mass matrix can then be 
expressed as 

M ac = {K p W)^H hc = K p YF ah E hc . (13) 

In the rest of this Section we calculate this mass matrix and analyze its eigenvalues for the 
three possible types of vacua: electroweak-symmetric, charge-breaking and neutral. 

3.2 Electroweak-symmetric vacuum 

The masses of the Higgs bosons in the electroweak-symmetric vacuum are determined only 
by the quadratic term of the potential and can be easily calculated in a straightforward way. 
The eight masses are grouped into two quartets with values ([T]) are 



m\ 2 = i i ( -nr n ) + (-m 22 



) ± ^ (m 2 u - m 2 22 ) 2 + 4\m 12 A . (14) 
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These masses squared are positive, if m\ x < 0, m| 2 < and m^m^ > |m 12 | 2 . However, we 
find it useful to work out this simple case in the reparametrization-covariant formalism just 
to illustrate how it works. 

The hessian H ab comes only from the M M r M term of the potential and is equal to — M^E^,. 
The mass matrix is then 

M ab = K p (-M p )(m p ) ah . (15) 

Matrices S's have a block-diagonal form, and therefore so does the mass matrix ffT5l) . It is 
built of two identical 4-by-4 blocks (A4 4 ) a b, with a,b = 1,2,3,4, whose form is still given by 
the same expression but now with 4-by-4 matrices S' 4 . In order to find its eigenvalues, let us 
calculate the trace of its successive powers: 

Tt[M 4 } = K p (-M M )Tr[S p S^] = -A(KM) , 

Tr[(_M 4 ) 2 ] = K p (-M p )K p/ (-M p/ )Ti[m^'^'} 

= 2(KM)K p M p r Ti[m p ] - K p K p ,M 2r Ti[m p '} = 8(KM) 2 - AK 2 M 2 , 

Ti[(M 4 ) n ] = -2(KM)Tr[(A<4) n_1 ] - K 2 M 2 Ti[{M 4 ) n - 2 ] . (16) 

These relations among the traces prove the mass matrix has only two independent eigenvalues. 
A simple analysis shows that there are two pairs of different eigenvalues, which are equal to 

m \ 2 = -(KM) ± a/ (KM) 2 - M 2 , (17) 

where we used K 2 = 1. This expression is reparametrization- invariant and can be calculated 
in any frame. In particular, in the original frame, where K p = (1, 0, 0, 0), we obtain 

m 2 12 = -M ± \M\ . (18) 

Using the definition of M M , one can immediately recover ( fl4l) . Eq. (fl8l) also shows that in 
order for the EW-symmetric extremum to be minimum, the four-vector M M must lies inside 
the backward lightcone. 



3.3 Charge-breaking vacuum 

Let us now find the mass matrix of the general 2HDM in the case of a charge-breaking vacuum. 
The hessian has the following form: 

H bc = 2k pv Yf hh ,Yy cd i£vLp cl . (19) 

All fields here must be understood as v.e.v.'s (y? a ), but to keep the notation simple, we will 
suppress the brackets. Thus, the mass matrix can be written as 

Ms = 2K p A li „Z»V(<p ® <p)Z u . (20) 

By construction, this is a 8-by-8 matrix. However, we know that it must have four flat 
directions corresponding to the Goldstone modes. We shall now get rid of these four flat 
directions by showing that there exists a 4-by-4 matrix A4 4 such that trace of any power of 
A4 8 is equal to the trace of the same power of A4 4 . 
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Indeed, consider the trace of A4g. Thanks to the properties of S's, we have 

Tt[M 8 ] = 2K p A^<pZ»Z^<p = 2K p A^(gWm" + gV-g^m p ) 
= 2K p m p (2A PP - TrA g pp ) = 2Tr [S ■ A] . 

Here, the matrix S ■ A is a symbolic form of the tensor S^ a A au = S^ a A aV) where 

= K v m^ + KV - (Km)g v ^ . (21) 

Note that the matrix S-A is defined in the euclidean space, and although it contains the tensors 
S 71 " and A au , they are contracted according to the usual rules of matrix multiplication. 
Consider now the trace of the square of 

Tr [(M 8 ) 2 ] = AKpA^KprA^ ■ ^E^V ■ <pZ u 'Z^<p . (22) 

Note that this expression does not factorize because A Mi/ and A p i u i couple the first and the 
second threads of S's. Consider one of these threads, e.g. 

This is a quadratic form in ip a ; therefore, only the afr-symmetric part of the product of S's 
survives. This effectively leads to the v <-> fx' symmetrization, and one can again apply ffl~2l 
to obtain 

The trace of the square of the mass matrix is then 

Tr [(M 8 ) 2 ] = AA^'A^rS^ = Tr [(2S • A) 2 ] . 

This calculation is easily generalizes to any power of the mass matrix: 

Tr l(M 8 ) n ] = Tr [(2S ■ A) n ] . (23) 

The fact that the trace of any power of A4g is equal to the trace of the same power of the 4- 
by-4 matrix 2S -A, means that there are four zero-modes in Aig and that all the four non-zero 
eigenvalues of coincide with the eigenvalues of 2S ■ A. Thus, the four eigenvalues of the 
matrix 2S ■ A gives the masses squared of the physical Higgs bosons in the charge-breaking 
vacuum. 

There is no simple way to calculate the masses themselves. However, the product of all 
four masses squared can be easily inferred from the above expression: 

Y[ m\ = det(2S ■ A) = 16 detS • detA . (24) 

i 

Both tensors here are written in the euclidean space. Determinant of euclidean A a p is the 
product of the eigenvalues^] of Minkowski A^: detA = A0A1A2A3. In order to calculate the 
other determinant, let us take a closer look at S^ v . The way it is defined, Eq. (I2T1) . allows 

2 Note a subtlety here: in a generic basis, the eigenvalues of the euclidean matrix A"^, which is not even 
symmetric, are different from the eigenvalues of the Minkowski tensor A**", i.e. Ao and Aj. However, the 
product of all the eigenvalues of these two matrices are equal. 
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us to immediately find its eigenvalues. Indeed, consider first a reduced version of this tensor, 
K^m u + K u m tl . In general, K^ and m fl are non-parallel four vectors, both lying strictly 
inside the forward lightcone. Within the subspace spanned by them, one can identify two 
eigenvectors of this reduced tensor, 

= 0, 




whose eigenvalues are (Km)±VK 2 m 2 . Note that e+ lies inside the forward lightcone, while 
lies outside it. In addition, there are two eigenvectors in the subspace orthogonal to K^ and m* 1 , 
with zero eigenvalues. Since adding a term proportional to g pv does not change the eigenvectors 
but just shifts all the eigenvalues by a common constant, we get the following result: S^ u is 
diagonalizable by an appropriate SO(l,3) transformation, and after diagonalization it take 
form: 



= diag(S , -S ll -S 2 , -S 3 ) , S = Vm 2 , S 1 = -Vm 2 , S 2 = S 3 = -(Km) . (25) 
Therefore, we obtain: 

\[m 2 = 16A (-Ai)(-A 2 )(-A 3 ) -m 2 (Km) 2 . (26) 

i 

As said above, a charge-breaking extremum exists, if m M lies inside the future lightcone, i.e. 
if m 2 > and (Km) > 0. It is also known that the charge-breaking extremum is a minimum 
if the tensor is positive-definite in the entire Minkowski space, i.e. if all its spacelike 
eigenvalues A 1j2j 3 are negativ^. Thus, all factors in fl26l) are positive. 

Another observation concerns cases when the potential has an explicit symmetry. Consider, 
for example, the lowest possible explicit symmetry, a Z2-symmetr;y@, which consists in reflection 
of, say, second axis. This explicit symmetry means that K 2 = 0, M 2 = 0, and that A 2fl = 
for /j ^ 2. It is known that the position of the charge-breaking minimum preserves all the 
discrete symmetries, so that m 2 is also zero. In this case one can immediately calculate the 
mass squared of the excitation that violates this symmetry: 

m 2 2 = 2(-A 2 )(Km) . (27) 



3.4 Neutral vacuum 

Let us now consider the neutral vacuum. The four-vector corresponding to a neutral 
vacuum must lie on the surface of the forward lightcone (again, we always refer to the v.e.v.'s, 
so that the brackets (• ■ •) are implicitly assumed). Therefore, the minimization procedure 
involves a Lagrange multiplier C, which brings up a new lightcone four-vector, defined as 
C M = A^ u r u — — ( ■ r M . This new four-vector gives rise to an additional term in the mass 
matrix: 

M 8 = 2K p K ilv t, p ^(<p ® ip)^ + K p Q^ p ^ . (28) 

3 We checked that these conditions can be also inferred from the positive-dehniteness of the mass matrix 
just derived. 

4 This symmetry is known in the literature as a generalized CP-symmetry. The "conventional Z2" corre- 
sponds, strictly speaking, to a (^2) 2 -symmetry of the potential, see details in [16] . 
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This matrix is again an 8-by-8 real symmetric matrix. However, one can easily split it into two 
4-by-4 matrices corresponding to the charged (the first four components of ip a ) and neutral 
(the last four components of (p a ) modes, which do not mix. 

Before we proceed, let us note that essentially this expression for the mass matrix of the 
most general 2HDM, but with a trivial kinetic part, was obtained in other works, [TT1 [T5l 18]. 
All these papers followed then the standard procedure: one switches to the basis where only 
the first doublet has non-zero v.e.v. (the Higgs basis), and then the entries of the mass matrix 
can then be written in a simple way via the parameters of the potential in this specific basis as 
well as v 2 . We show in this subsection that the basis- invariant features of the mass matrix can 
be written in an SO(l, 3)-covariant way, without referring to any specific basis. The power of 
the covariant expression is that it can be analyzed in any desired basis, e.g. in the A M „-diagonal 
basis. We checked that in the canonical basis, our result reproduce those of [151 IB] • 

Consider first the charged excitations. Their masses arise solely from the last term in fl28|) : 

Mt = K p (^ p ^ , (29) 

where S's are now 4-by-4 matrices. By explicit calculations and using the fact that C 2 = 0, 
we get: 

TrMt = 4(K{) , Tr[(M? f) = 8(K() 2 , Tr[(Mt) n ] = 2[2(KC)] n • (30) 

It means that this matrix has only two non-zero eigenvalues, which are identical and equal to 

m 2 H± = 2{K() . (31) 

This implies, in particular, that in order for the extremum to be minimum, £ must lie on the 
surface of the forward, not backward lightcone. 

For the neutral modes one has the same expression as in ( l28i) . but with 4-by-4 matrices 

E M : 

Ml = 2K p A^W( v <g> ^E" + K^TPW . (32) 
Let us calculate the trace of the mass matrix of the neutral Higgs bosons: 

TrMt = 1K p k^Y?Y, p Y»Lp + 4(K() = AK pu K»r v - 2TrA (Kr) + 4{K() 

= 2(4A^ - TrA 9ilv )K»r v - A(KM) . (33) 

We expect that among the four neutral modes there will be one goldstone, which makes the 
determinant of M%' zero. To check it explicitly, we first factor out the matrix K p Y, p and 
check by a direct calculation that its determinant is equal to (K^K^) 2 = 1. The remaining 
determinant 

det [2A^E»((p®ip)E u + (^] 

is equal to zero, which can be best seen in the Higgs basis, where the second row and the 
second column have only zeros. In the generic basis, the goldstone mode is Wi = (n o )jj0j, 
where the matrix n° is the generator of the SO (2) rotations between the real and imaginary 
parts, see Appendix. 
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3.5 The extra symmetry of the neutral modes 

The appearance of the tensor 4A MZ , — Tr A g pv in (1331) is not accidental, but reflects an extra 
symmetry of the neutral mass matrix. If we consider the neutral vacuum and if we analyze 
only neutral excitations, we always stay on the surface of the lightcone: we consider only r M 
such that g^r^r 1 ' = 0. This means that if we shift the tensor A^ in the potential as 

A Mi , -> A^ + CV (34) 

with an arbitrary C, the purely neutral contribution to the potential does not change, and 
neither does the neutral mass matrix. The tensor 4A M1/ — Tr A g^ v is precisely the combination 
that is invariant under such a shift. In terms of the original parametrization of the quartic 
potential ([I]), this symmetry means that the neutral Higgs boson masses do not depend on 
the value of TrA = A3 — A4. 

One can make use of this extra symmetry to simplify the neutral Higgs boson mass matrix. 
First, note that the neutral mass matrix ( 1321) is invariant under the transformation ( 1341) thanks 
to the following relation: 

2g^{ip <g> <p)Z v + = . (35) 

Let us recall now that is proportional to r M : = £ ■ r M , where ( is the Lagrange multiplier 
of the minimization problem. Then, we can group the two terms in (1321) together: 

Ml' = 2K p k^ p ^{ V ®v)Y, u where A^ u = - (g^ . (36) 

It is remarkable that the new tensor A M „ is itself invariant under (1541) as this shift is accom- 
panied by C — > ( + C: 

Cr^ = C = A^rv - M M -> (A^ + Cg^)r„ - M M = ^ + Cr^ 1 = (C + C)r» . (37) 

With this expression in hand, we can again use the trick from the analysis of the charge- 
breaking vacuum and state that all the neutral boson masses are given by the eigenvalues of 
the following matrix written in a manifestly covariant form: 

Ml = 2S ■ A , where S" v = K"r u + K v r p - {Kr)5» . (38) 

Therefore, one can immediately write the trace of any power of the mass matrix: 

Tr[{Mt) k ] = 2 fc ^VA% 2 • • • S'^A'V , (39) 

and calculate the determinant of S^ v using (125]) : 

det£ = -r 2 [Krf = , (40) 

which proves the existence of a goldstone mode in a basis-invariant fashion. 

4 Discussion and conclusions 

The principal result of this paper is a demonstration that the mass spectrum of the general 
2HDM can be studied in a reparametrization-invariant way within the Minkowski-space for- 
malism of [T5l [16] . This means that the scalar propagators can be now written explicitly and 
can be used, for example, to improve the thermal one-loop calculations of [20] . 
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Another interesting issue that one can now address is to understand to what extent the 
perturbativity/tree-level unitarity bounds on the Higgs potential restrict the values of the 
Higgs boson masses. In the Standard Model, there is a strong correlation between the value 
of the quartic coupling constant A and the Higgs boson mass. Therefore, an upper limit 
on A implies a corresponding upper limit on Mg. In the 2HDM, due to a large number of 
free parameters, the situation is more complicated, see [221 ffl EH E31 EI] • It was noted that 
in certain cases masses of some of the Higgs bosons can be very high without violating the 
tree-level unitarity conditions. With an explicit expression for the trace of the mass matrix, 
one could now attack this problem in the most general case within the Minkowski-space 
technique. The only piece still missing is a reparametrization-covariant expression for the 
tree-level unitarity constraints. 

In conclusion, we showed that the Minkowski-space approach to the most general 2HDM 
can also be used to analyze the mass spectrum of the physical Higgs bosons. We calculated 
the traces of the powers of the mass matrix and its determinant for all types of vacuum that 
can exist in 2HDM. These results can now be used to get even more insight into the properties 
of the general 2HDM. 

We are thankful to J.-R. Cudell for helpful discussions. This work was supported by the 
Belgian Fund F.R.S.-FNRS via the contracts of Charge de recherches (I.P.I.) and of Aspirant 
(A.D.). The work of I.P.I, was in part supported by grants RFBR 08-02-00334-a and NSh- 
1027.2008.2 



A Algebra of matrices and II 



The four- vector of matrices E^ 1 is introduced via Eq. (Q. The full 8-by-8 matrices S M have 
block-diagonal form and are built from two identical 4-by-4 matrices, which we also denote by 
the same letter S's and whose properties we describe here. 

S° is just the unit matrix, while the explicit expressions of are: 
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o \ 
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(41) 



These matrices satisfy the Clifford algebra condition: 

{p,TP} = 25 ij I 4 



(42) 



The set of S's is not closed under taking commutators. Instead, they can be expressed via 
real antisymmetric matrices IV: 



( 



it = n°s 



where n 



o 



V 



o 
o 



i 















1 

0/ 



(43) 



The matrix n° is the generator of the simultaneous 150(2) rotations between the real and 
imaginary parts the two doublets; it commutes with all £ l and its square is equal to —1. The 
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set of matrices S l and II* now forms the algebra: 

[£*, £ J '] = 2e ijk U k , [£*, n J '] = -2e^' fc £ fe , [IT, H" J '] = -2e^ fc n fc . (44) 

Note that II 1 do form a closed algebra. 

The algebra of S* and II* is isomorphic to the usual Poincare algebra of the generators of 
boosts and rotations. Using this, we can introduce matrices 

X l ± = \ (±E j - iW) , (45) 

which satisfy the following commutation laws: 

[X t ± ,xL] = uP k X$ i1 [X t ± ,Xi F ]=0. (46) 

Finally, we note that any four-vector a M can be associated with a real symmetric matrix 
A = a M S M , which has the following properties: 

detA = (aX) 2 , A ~ l = , with s " = • (47) 
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